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The magnetic moment of tau neutrinos in the MeV mass range may be 
large enough to modify the cosmological freeze-out calculation and determine 
the tau-neutrino relic density. In this paper we revisit such a possibility. We 
calculate the evolution and freeze-out of the tau neutrino number density as 
a function of its mass and magnetic moment. We then determine its relic 
density, then calculate its effect upon primordial nucleosynthesis including 
previously neglected effects. 
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Present experimental bounds to the electromagnetic properties of the tau neutrino 
are several orders of magnitude less stringent than the bounds to the corresponding 
properties for electron and muon neutrinos. For instance, while the upper limits to 
the diagonal magnetic moments of the electron and muon neutrinos are is a Bohr 
magneton) 1.1 x 10~ 9 /i-B and 7.4 x 10~ 10 /xb respectively JT|, the experimental upper limit 
to the diagonal magnetic moment of the tau neutrino is \l Vt < 5.4 x lO~ 7 [i B 0. More 
stringent bounds on neutrino magnetic moments of order 10~ 10 to 10~ 12 are available from 
astrophysical constraints ||, mainly from the cooling of stars and from the study of SN 
1987A. However, these bounds apply only if mass of the neutrino species does not exceed 
the stellar temperatures relevant for neutrino production. Furthermore, astrophysical 
constraints are model dependent since they assume that the outgoing wrong helicity 
neutrinos are completely sterile [[|. 

Big-bang nucleosynthesis (BBN) is a precious tool that has been employed to con- 
strain many neutrino properties so it is no surprise that BBN can be used to bound 
neutrino magnetic moments. In 1981, Morgan || showed that the "sterile" right-handed 
degree of freedom of Dirac neutrinos^ can be populated through processes like eui — > evR 
and e"e + — > v R v L , mediated through a virtual photon coupled to the neutrino through 
its magnetic moment. The degree to which the right-handed neutrino is populated de- 
pends upon the strength of the above interactions, which in turn is proportional to the 
magnitude of the neutrino magnetic moment. Thus, endowing a neutrino species with 
a magnetic moment potentially leads to an unacceptable doubling of the contribution 
of that species to the energy density, jeopardizing BBN's successful predictions. If the 
neutrino species is relativistic at freeze out, one must require that right-handed neutrinos 
decouple before the QCD phase transition so that their number density is diluted by the 

1 If the neutrino has a magnetic moment it must be a Dirac fermion (we do not consider transitional 
magnetic moments in this paper). 
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huge entropy shift associated with the transition. Such a requirement translates into 
an upper limit to the neutrino magnetic moment: fi Ur < 1 to 2xl0~ 11 /i£. However, in 
deriving this limit Morgan did not consider the possibility of non-zero neutrino masses. 
There are two main effects to be considered if one allows a non-zero mass for the neu- 
trino species in question. First, the neutrino may not be relativistic at freeze out, and 
its number density must be calculated by solving the Boltzmann equation. Furthermore, 
the scaling of the neutrino energy density with temperature depends upon the neutrino 
mass. Therefore, Morgan's useful limit applies only if the neutrinos are ultra-relativistic 
around freeze out and BBN, that is if m VT < O.lMeV. 

The tau neutrino could be heavier than 0.1 MeV;0 in fact, the present upper bound 
to the tau-neutrino mass is m„ T < 24 MeV [0. Although somewhat model dependent, 
more stringent bounds on the neutrino masses can be determined from cosmological con- 
siderations. In particular, if the relic heavy-neutrino energy density today is sufficiently 
large, the predicted age of the universe will be less than observed. If the neutrino is 
stable and it is nonrelativistic today, the age limit (Qh 2 < 1) constrains the mass of any 
stable neutrino species to be less than the Cowsik-McClelland limit, m v < 91.5 eV 
Of course if the neutrino is unstable, the Cowsik-McClelland limit can be evaded ||. 
But even in this case there are lifetime-dependent limits to the neutrino mass. If the 
heavy-neutrino lifetime is longer than a second or so, it can give an additional contri- 
bution to the energy density during nucleosynthesis and spoil the successful predictions 
of standard calculations. Using these kind of considerations, BBN constraints to the 
tau-neutrino mass excludes the range 0.3 < m u < 25 MeV if it is a Dirac fermion, and 
the range 0.5 < m v < 25 MeV if it is a Majorana fermion [K|. (It was assumed in the 
above BBN analysis that the neutrino eventually decays after BBN; if it decays after 



2 We assume here that the mass of the muon neutrino is less than 0.1 MeV. 
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decoupling but before or during BBN, the situation is more complicated pT| .) 

The Cowsik-McClelland limit and the nucleosynthesis considerations might be mod- 
ified if one introduces new interactions that changes the neutrino annihilation cross sec- 
tion. This is the case if the neutrino has a large diagonal magnetic moment, because a 
large magnetic moment would increase v—v annihilation (creation) into (by) e , keeping 
the neutrinos in equilibrium below the canonical (including only weak processes) neu- 
trino decoupling temperature of about an MeV. If the neutrino mass is sufficiently small 
(much less than an electron mass) and remains coupled to electrons while the electrons 
annihilate, the neutrino number density will be increased because part of the electron's 
entropy will be shared with the neutrinos. However, if the neutrino mass is not much 
less m e and it remains in equilibrium through magnetic-moment mediated interactions, 
its energy density will be Boltzmann suppressed before decoupling, weakening the BBN 
constraints. 

In this letter we study how the interplay between the neutrino mass and magnetic 
moment modifies the cosmological constraints to tau-neutrino properties from the age of 
the universe and BBN. Besides the mere extension of the upper limit on //„ T to larger 
neutrino masses, the main purpose of our letter is to give a final answer to the intriguing 
possibility that tau neutrinos with a large magnetic moment could form cold dark matter. 

Giudice [|l2j first observed that if tau neutrinos are stable, have a mass in the range 



m„ T ~ 1 to 10 MeV, and are endowed with a magnetic moment of \i Vt ~ 10 _6 /i£, they 
would stay in equilibrium through their magnetic-moment interactions and would decou- 
ple when they are nonrelativistic. If the magnetic moment is large enough, their final 
abundance might give rise to a universe with VL v h 2 ~ 1. Although the latest experimental 
upper limit on fi Ur |fj] seems marginally at odds with Giudice's scenario, it is worthwhile 
to investigate this hypothesis further . 

Giudice made use of the fact that Morgan's conclusions about doubling the effective 



tau-neutrino number density by populating the right-handed component can not be ap- 
plied directly to MeV-mass neutrinos. This is because their energy density is Boltzmann 
suppressed at freeze-out and during BBN. Thus, even including the right-handed compo- 
nents, tau neutrinos will not contribute so much to the energy density as to spoil BBN. 
We show that while this is approximately true, it is not exactly true. BBN is such a 
sensitive probe of the expansion rate of the universe at the temperatures of interest that 
even a small contribution to energy density is important. Therefore the contribution of 
the right-handed neutrino to BBN requires a careful treatment. We report the results of 
such an investigation in this communication. 

There are three effects that must be carefully accounted for: 1) After a massive 
neutrino species decouples and becomes nonrelativistic its energy density grows rela- 
tive to the energy density of a massless neutrino species fl4" |. Although one must solve 



the Boltzmann equation to compute the energy density of the heavy neutrino (see be- 
low), it is possible to estimate this effect by observing that p v {m v ^ 0)/p v (m v = 0) ~ 
(m u /3.\5T u )r oc t 1/2 when T ~ m V} where r is the ratio of the number density of massive 
neutrinos to massless neutrinos after freeze-out. 2) A neutrino species with a mass in 
the MeV range and with a magnetic moment close to the present experimental limit de- 
couples when it is semi-relativistic. Neither the relativistic nor the nonrelativistic cross 
section can be used, and a general treatment of the thermal-averaged annihilation cross 
section used in the Boltzmann equation for the neutrino abundance is required. (Giudice 
performed his analysis in the extreme nonrelativistic limit.) 3) Plasma effects must, at 
least a priori, be considered. Not only do thermal corrections to the amplitudes of the 
main processes involved in BBN have to be included ||15|| , but more importantly, the mass 
corrections due to the electromagnetic coupling of the particles to the relativistic plasma 
must be accounted for. For example, the photon in the thermal bath becomes a plas- 
mon and acquires an effective mass p6|. The plasmon mass has a double effect. In first 



place, it affects the electromagnetic channel of the neutrino annihilation cross section 
Although this is a second-order effect some resonance might enhance it dramatically []IT 



Secondly, a plasmon mass gives rise to new processes that are kinematically forbidden in 
the vacuum. In our case the most relevant of these processes is the decay plasmon — > vv, 
having a rate [18 

r vh T ( 2 A 2 \ 3 / 2 K±{x P ) 

where xp = rrip/T with mp the temperature-dependent plasmon mass, up ~ 0.1T is 
the plasma frequency, and Ki(x) are the modified Bessel functions of order i. However, 
the threshold up > 2m v reduces the importance of the process plasmon — > vv during 
BBN for Me V- mass neutrinos. Analogously, since 2m l/ > wp, screening effects induced 
by the plasma on the photon propagator turn out to be negligible. The relative unimpor- 
tance of these considerations were verified by directly including them in our numerical 
calculations. 

The cross section for the electromagnetic channel of the process vv — > e~e + is 

= ^ (r^lj;) (l + 8 ^ + 2 ^ + 16 ^) , (2) 

where \fs > 2m e is the total center-of-mass energy. 

The weak contribution to the annihilation process of Eq. (Q) can be neglected if the 
neutrino magnetic moment is larger than 10 _10 (m jy /l MeV) /zb- We will work within the 
limits of this assumption.^ 

Because both helicity eigenstates of the neutrino are symmetric with respect to elec- 
tromagnetic interactions, we do not differentiate between them in our calculations. For 

3 If the magnetic moment is larger than 10 _lc \ra v / '1 MeV)/is, then neutrino annihilation will occur 
predominantly through photon exchange, rather than Z exchange. For m„ < 100 keV, considerations 
of stellar energy loss by neutrino pair emission limits the magnetic moment to be greater than about 
10~ n /XB. Thus, we will consider neutrinos more massive than 100 keV. 
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this reason the processes ei>L{R) i — ► ^r{l) changes neither the total, nor the relative vl 
vs. vr abundances.^] 

The Boltzmann equation for the abundance of the heavy neutrino is [03 

dY ( 7i \ V 2 gl /2 m u m P i / 2 2 x 

where x = m u /T, Y = n v /s is the ratio of the v T number density to the total entropy 

— 1/2 

density of the universe, vyy^ is the M0ller invariant flux factor, and m Pi = G N is the 
Planck mass. The parameter g* is defined as 

g ;/2 = ^( 1+ i T dh tf/ m \ 



1/2 
9eff 



3 h eff (T) dT 

where the effective number of degrees of freedom for the energy density, g e ff(T), and for 
the entropy density, h e ff(T), are defined as 

7T 2 27T 2 

P = 9eff(T)-T^ s = h eff (T)—T* . (5) 



Following [B5|, the thermal averaged cross section times the M0ller velocity is 

1 f°° 

{aVM ^ = SmiTKKx) L v a{3) (S ~ 4 -')^i(v^/T) ds. (6) 
We have used the Maxwell-Boltzmann distribution to compute the thermal-averaged 
cross section (for a detailed review of computations in this approximation see e.g., Ref. 
|20|). Although normally this is a very good approximation only for temperatures 



T <3m u , we have checked that at the freeze-out temperature (the only temperature 
around which Eq. (P) plays a relevant role) the approximation is adequate. 

The neutrino decoupling temperature Tp is here defined by the condition Y(Tp) — 
Yeq(T f ) = 1.5Y EQ , where 

_ njjQ _ 45 I v {x) 
S 7T 4 h eff (T) 



4 The vl(r) helicity eigenstates should not be confused with the chirality eigenstates. Since we ignore 
weak interactions, chirality does not play a role in our analysis. 
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Figure 1: Tau neutrino abundance vs. the parameter x = m VT /T is represented for dif- 
ferent value of \l Vt : from below, the three different curves refer to \l v = 10~ 6 /ib, 10~ 7 //b, 
and 10~ 8 /j,b- Here we have chosen m u = 1 MeV. The logarithms are base 10. 



We numerically solved Eq. (§) to compute the tau-neutrino abundance as function of 
T for fixed values of ra Vr and fi„ T . Since the freeze-out temperature increases with \x Vt , it 
is natural to expect that the final tau-neutrino abundance is suppressed as the magnetic 
moment increases. This is clearly visible in Fig. 1. Assuming the tau neutrinos to be 
stable we can easily check their effect on the dynamics of the universe for several values 
of m Vr and \l Vt . In particular, we first consider the contribution to the present energy 
density of the universe due to massive tau neutrinos: 



with 





tt UT h 2 = 




(9) 



1.054 MeVcm- 3 ' 
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Figure 2: The effective number of degrees of freedom in the energy density is shown as 
function of temperature for two chosen value of fi v and m v = 1 MeV. The upper solid 
curve is for /j, Ut = lO~ 8 fi B ', the lower solid curve is for /jl Vt = 10~ 6 //b- The reader can 
compare our result with the result obtained for 3.4 standard massless neutrinos shown 
by the dashed line. The logarithm is base 10. 

where indicates quantities evaluated at the present time. Requiring VL Vt Ii 2 < 1, we can 
verify which region of the parameter space /i Ut versus m Vr is compatible with the age 
constraint. 

Of course if the tau neutrino is unstable, the cosmological age constraint discussed 
above does not apply. However we can still use BBN to limit the properties of the tau 
neutrino provided that the lifetime, T„ r , is greater than about a second. 

To evaluate the impact of the massive tau neutrino with a large magnetic moment 
on BBN we must know how it modifies the effective number of degrees of freedom of the 
energy density, g e ff(T), for 0.1<T<10 MeV, since the light element relic abundances 
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depend critically on the expansion rate of the universe during BBN, which in turn is 
parameterized by g e ff 0: 

H(T) = l.mg l J f ){T)^- . (10) 
The tau neutrino contribution to g e ff is given by 

g„AT)=P„ T (T)(^)T-\ (11) 

where p Vr = sY Ut yj (3.15T Ut ) 2 + , and T Ut is computed by imposing entropy conserva- 
tion. Fig. 2 clearly demonstrates that g e ff grows as the decoupled tau neutrinos become 
nonrelativistic. This effect becomes less pronounced as the magnetic moment is increased. 
Of course this is simply because increasing \x Vt decreases T F , leading to a tau-neutrino 
energy density more effectively Boltzmann suppressed before freeze out. For this reason, 
values of the tau-neutrino magnetic moment larger than 10 -8 \xb are not expected to have 
a large effect on BBN if m Vr > 0.1 MeV. 

To check this in detail and in order to be able to evaluate the effects of the neutrino 
mass and magnetic moment on light element production, we incorporated our results for 
the abundance as a function of temperature into the standard nucleosynthesis code pTf . 
In Fig. 3 our predictions for the relic 4 He abundance as a function of the tau-neutrino mass 
are shown for two values of the magnetic moment. As expected, the predicted abundance 
Yp is suppressed with increasing \i Vt . Increasing the mass above a few MeV increases Yp, 
since the tau neutrinos then become nonrelativistic earlier. For small masses, Yp grows 
with decreasing m VT . This is due both to the less effective Boltzmann suppression and 
to the entropy transfer from e ± annihilation to the tau neutrinos. 

In order to discriminate which region of the m u versus p u parameter space is com- 
patible with observations, we require that the predicted light element abundances do 
not exceed the observational limits [|2j: Y P < 0.24; (D + 3 He)/H < 1.1 x 10~ 4 ; and 
7 Li/H < 1.7 x 10~ 10 . Since the baryon-to-photon ratio rj is a free parameter, for ev- 

9 



0.5 1 



5 10 

m (MeV) 



Figure 3: The predicted 4 He relic abundance is represented as function of the tau-neutrino 
mass for two values of ji v . The upper curve is for /i Ut = 10~ 8 {i B , while the lower one is 
for h Vt = 1CT 6 fiB- The dashed line corresponds to the observational upper limit. 

ery chosen pair of m Vr and {i Ut we fix it at the minimum value compatible with the 
(D + 3 He)/H upper limit. Then we check if the predicted 4 He relic abundance is consis- 
tent with the upper limit of 0.24. The 7 Li constraints turn out to be always less stringent 
than the limits coming from 4 He. 

Our results are summarized in Fig. 4. As the reader can observe, the age-based 
constraints are much more stringent than BBN constraints if the tau neutrino is stable. 
In this case the border between the allowed and forbidden regions in the m VT versus. ji Ur 
parameter space from the age constraint almost coincides with the experimental limit 
line. This is a remarkable coincidence. In fact, the age limits are lower limits to (j, Ut , 
whereas the experimental limits are upper limits. This means that nearly the entire 
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Figure 4: Exclusion plot in the tau neutrino magnetic moment — mass parameter space. 
The solid line provides the lower limit to /i„ T coming from the requirement Vth 2 < 1. 
The dashed line provides the corresponding limit from BBN considerations. The dotted- 
dashed line represent the experimental upper limit. 

parameter space for O.lMeV < m. Ur and 10 _10 (m^ T /lMeV) < /i„ T (the very range to which 
our consideration apply) is excluded by our considerations. To be precise, a very small 
region between the experimental and the age-based constraints remains open. Even this 
region would be closed using a slightly larger value of h as recent observations suggest. 

As a consequence, Giudice's hypothesis is definitely ruled out. Furthermore, our 
results improve the upper limit on the tau-neutrino magnetic moment by several orders of 
magnitude in the mass range we considered. We have checked that plasma-physics effects 
are subdominant. We have to stress that this limit is valid only if the tau neutrino is 
stable, as indeed Giudice assumed. Stability can be achieved by imposing some additional 
symmetries, e.g., individual lepton-number conservation. Of course, in any case some new 
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physics beyond the standard model must be introduced in order to have such a large value 
of \x Vt . Furthermore, a tau neutrino with mass larger than 1.1 MeV decaying according to 
the a minimally extended standard model via the channel v T — > u e e + e~ is incompatible 
with BBN. In fact, since experimental data constrains the v T lifetime to be 1 s< t Vt < 10 
s, if itl Vt > 1.1 MeV, electrons and positrons produced from this decay would induce the 



photodestruction of light elements [23 



If the tau neutrino is unstable but the lifetime exceeds one second, a band of magnetic 
moment values, roughly 10~ 8 /ib < ^v T < 10~ 6 /i_B, remains compatible with experimental 
and cosmological bounds. This confirms the result of Ref. j24] and extends it to a wider 
tau-neutrino mass range. It is understood that in this case the tau neutrino has to decay 
in some non-standard way in order the decay products do not affect dramatically the 
light element relic abundances. 
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